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Abstract. When astronomers study the dark matter halos of spiral galaxies, they normally assume that the disk mass-to-light 
ratio is constant. We describe a method of analyzing the kinematics of planetary nebulae (PNe) in nearby face-on spiral 
galaxies to test this assumption. Since the restoring force for stellar motions perpendicular to the galactic disk is proportional 
to the disk mass surface density, measurements of the vertical velocity dispersion can be used to produce an independent 
measure of the total amount of matter in the disk. Our steps are: (1) to identify a population of PNe by imaging the host spiral 
in several filters, and (2) to isolate the vertical velocity dispersion from spectroscopic observations of the PNe. Our first results 
for the PNe of M33 indicate that the mass-to-light ratio of the galaxy's disk actually increases by more than a factor of 5 over 
the inner 6 disk scale lengths. We have begun similar studies of the PNe in five more face-on galaxies: M83, M101, M94, 
NGC 6946, and M74. These data will also produce additional science such as galaxy distances and constraints on the disk 
transparency. 



INTRODUCTION AND SCIENTIFIC 
MOTIVATION 

Dark matter, a mysterious topic, has been under intense 
study for many years. One way to shed some light on the 
subject is to examine galactic halos. In the case of spiral 
galaxies with weak bulges, one typically measures the to- 
tal galactic mass via the system's rotation curve (Faber & 
Gallagher 1979; Ashman 1992; Combes 2002), assumes 
that the mass near the center of the galaxy is entirely 
baryonic, and subtracts off the contribution of the disk by 
assuming a constant disk mass-to-light ratio (Kent 1986, 
Palunas & Williams 2000; Sofue et al. 2003). (This is 
the "maximal disk" method.) However, while absorption 
line studies (Bottema 1993; Bottema, van der Kruit, & 
Freeman 1987; van der Kruit et al. 2001; Gerssen, Kui- 
jken, & Merrifield 1997, 2000) have indicated that the 
disk mass-to-light ratio is constant in the inner regions 
of spiral galaxies, there formerly were no results about 
the outer regions where the influence of dark matter is 
greatest. The extreme difficulty in separating the mass of 
a galaxy's visible disk from that of its dark halo limits our 
understanding of almost every facet of galaxy formation. 
An independent method of determining the disk mass is 
needed to break the disk-halo degeneracy. 

OUR METHOD: USING PLANETARY 
NEBULAE TO STUDY DISK MASS 

To determine a disk's mass-to-light ratio, one can study 
the vertical motions of old disk stars. Since the restor- 
ing force for stellar motions perpendicular to a galactic 



disk is proportional to the disk mass surface density (£), 
measurements of the vertical velocity dispersion (ct z ) im- 
mediately yield an independent measure of the amount of 
matter in a disk. From the isothermal disk approximation, 
old disk stars oscillate in z according to 

aj{R) = nGZ(R)z , (1) 

where zq is the scale height of the stars (Binney & 
Tremaine 1987). Since studies of edge-on spirals demon- 
strate that zo is constant with radius (van der Kruit & 
Searle 1981; Bizyaev & Mitronova 2002), this parameter 
can be fixed at some appropriate value. Then, by observ- 
ing a face-on galaxy, we can see if the matter scale length 
does indeed decline in a manner similar to the light. 

The traditional method of determining ct, via absorp- 
tion line spectroscopy is very challenging both observa- 
tionally and in terms of data analysis. Moreover, because 
the technique is limited by surface brightness, its effec- 
tiveness is restricted to a galaxy's inner regions (Bottema 
1993; Gerssen, Kuijken, & Merrifield 1997, 2000). An 
alternate method to measure a z is to use planetary neb- 
ulae (PNe) as kinematic test particles. PNe are relatively 
numerous, easy to detect in a galaxy's outer regions 
where dark matter is most important, and come from 
a progenitor population of low and intermediate mass 
stars. They are therefore representative of the old stellar 
disk of a galaxy. In addition, since PNe are strong emis- 
sion line sources, virtually every object that can be found 
photometrically can be observed spectroscopically. With 
a medium (R ~ 5000) resolution instrument, radial veloc- 
ity measurements to a precision of ~ 2 km s~ 1 are obtain- 
able without much difficulty. This makes PNe ideal test 
particles for probing the disk mass of face-on spirals. 



FINDING THE PNE 

In order to study the disk mass-to-light ratio, we must 
first identify a suitably large population of PNe. The spi- 
ral galaxy under study must be imaged with a 4-m class 
telescope in four filters: two narrow on-band filters cen- 
tered at [O III] A5007 and Ha (in the rest-frame of the 
galaxy), and two wider off-band filters, such as V and 
R. Since PNe have virtually no continuum, they can be 
found by blinking the on-band images against their off- 
band counterparts. True PN candidates are: (1) consistent 
with point sources, (2) detected in A5007 but invisible in 
V and R, and (3) invisible or weak in Ha. These criteria 
exclude most or all H II regions and supernova remnants. 
At this stage, we can determine the distance to the galaxy 
via the Planetary Nebula Luminosity Function (Ciardullo 
et al. 2002), and attempt to constrain the transparency of 
the disk via PN number counts. Thus far, as part of our 
program, we have discovered 152 PNe in the Triangu- 
lum Galaxy, M33 (Ciardullo et al. 2004; hereafter C04), 
65 PNe in the Northern Pinwheel Galaxy, M101 (Feld- 
meier, Ciardullo, & Jacoby 1996), and ~ 200 PNe in the 
Southern Pinwheel Galaxy, M83 (Herrmann, Ciardullo, 
& Vinciguerra 2005). 

PNE KINEMATICS 

Once we have identified a large population of PNe and 
determined accurate positions, the next step is to obtain a 
high precision radial velocity for each object. This can be 
done with fiber-coupled spectrographs, such as the Hydra 
instruments on the WIYN and Blanco 4-m telescopes, 
and the Medium Resolution Spectrograph on the Hobby- 
Eberly Telescope. In order to minimize systematic errors, 
it is useful to target the PNe multiple times and to pay 
special attention to the wavelength calibration. Ideally, 
the velocity uncertainties should be < 5 km s _1 . 

The top frame of Fig. 1 shows the preliminary radial 
velocities of 203 PNe in M83 after correcting for the 
systemic and barycentric velocities. PNe with the largest 
velocities away from us are clustered in the lower left of 
the frame, while PNe with the largest velocities toward 
us are in the upper right. Clearly, we are detecting the 
rotation of the galaxy. For a flat axisymmetric system, 
the radial velocity, v rad is given by 

v rad = V0cos0sini + v/jsin0sini + v z cosi, (2) 

where vq, vr and v z are the azimuthal, radial and verti- 
cal components, ;' is the inclination, and is the angle 
from the principal axis in the galaxy plane. The rotation 
component is eliminated by subtracting out v ro tCos sin; 
where v rot is taken from studies of H I gas (Tilanus & 
Allen 1993). The resulting residual velocities are random 
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FIGURE 1. PN radial velocities in M83. From fastest away 
to fastest toward the symbols are o, □, o, +, x , #. PNe not ob- 
served are indicated with a dash. The top panel shows the radial 
velocities as observed; the bottom illustrates the velocities after 
the removal of galactic rotation. 



with position in the galaxy, indicating that the rotation 
has been taken out. (See the bottom frame of Fig. 1.) 

If the line-of-sight velocity measurements are dom- 
inated by the vertical component of the velocity ellip- 
soid, then the dispersion in the residual velocities, G ra d, 
should be related to the disk mass surface density, E, via 
Eqn. (1). For a constant mass-to-light ratio disk, this im- 
plies that o rad should fall off exponentially with a scale- 
length twice that of the galaxy's light. Figure 2 shows the 
line-of-sight residual velocity dispersion for the PNe in 
M33 and M83. For M33, the dispersion clearly does not 
follow the light. This is not entirely unexpected, since the 
galaxy's ~ 56° inclination ensures that all three compo- 
nents of the velocity ellipsoid contribute to a rad . Our pre- 
liminary results for M83 are similar. Thus, in both cases, 
we need to de-couple a z from the other components of 
stellar motion. 
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FIGURE 2. Dispersion in the residual PN radial velocities for M33 and M83 after galactic rotation has been removed. The curves 
show the exponential decay expected for constant mass-to-light ratio disks. The error bars in y represent 90% confidence intervals, 
while those in x indicate the bin size. 



The observed line-of-sight velocity dispersion of a 
galaxy can be written in terms of the azimuthal (<j>), radial 
(R), and perpendicular (z) velocity dispersions via 

°md = cos 2 sin 2 ;' + a| sin 2 ^) sin 2 ;' + a? co s 2 ;' + a 2 ea , , 

(3) 

where G meas is the measurement uncertainty. (Note that 
to extract G z from G ra d, G mea s must be kept to a mini- 
mum.) Since there are three unknowns but only one equa- 
tion, external constraints are needed. One such constraint 
is the epicyclic approximation, which allows us to re- 
move G§ from the equation by writing it in terms of Gr 
and the radial gradient of the circular velocity (Binney 
& Tremaine 1987). Two others are the Toomre (1964) 
criterion, which requires that the disk be stable against 
axisymmetric perturbations, and the firehose instability, 
which forces us to consider only those disks that are sta- 
ble against buckling (Toomre 1966; Merritt & Sellwood 
1994). Finally, for non-barred galaxies, we can impose 
the Morosov (1980, 1981a, 1981b) criterion, which re- 
quires that a disk be stable against the formation of a bar. 

Each of these constraints eliminates some combina- 
tion of G z and Gr from consideration. We can then use 
a maximum-likelihood analysis to determine which of 
the remaining combinations of the two variables are most 
probable. If enough PNe are observed, we can also deter- 
mine the system's asymmetric drift, and place a further 
constraint on the solution (Binney & Tremaine 1987). 



RESULTS FOR M33 

At present, M33 is the only galaxy for which we have a 
complete analysis. The most likely values of G z and Gr 
are shown in Fig. 3. The data demonstrate why the line- 
of-sight velocity dispersion of M33 varies so little with 
radius. Near the center of the galaxy, the increase in G z is 
negated by a turnover in the radial velocity dispersion. (If 



this did not happen, then Gr in the central ~ 1 kpc would 
be greater than v ro t, and the galaxy would have a bulge.) 
At larger radii, both Gr and G z decline exponentially, 
with the scale length of Gr being ^25% larger than 
that of G z . The derived values of the dispersion ratio, 
as well as its radial gradient, are in excellent agreement 
with models of disk heating (Villumsen 1985; Jenkins & 
Binney 1990; Carlberg 1987). 

A more surprising result is that the mass scale length 
inferred from G z is more than twice that expected from 
observations of the galaxy's infrared light. Near the cen- 
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FIGURE 3. The vertical and radial dispersions of M33 PNe 
derived from a maximum-likelihood decomposition of G mc i- 
The non-Gaussian uncertainties are illustrated by the width of 
the lines. The dotted curves show the best-fit exponentials. 
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FIGURE 4. The mass-to-light ratio of M33's disk as a func- 
tion of galactic radius. The error bars do not include systematic 
errors. 



of five additional galaxies which are more face-on than 
M33, including M83, M101, M94, NGC 6946, and M74. 
These data should tell us whether the spiral disks are 
truly super-maximal, as the M33 data suggests, or if Or 
and a§ are conspiring against us. 
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ter of the galaxy, the disk mass-to-light ratio of M/Ly ~ 
0.3 agrees with that derived from the galaxy's rotation 
curve. However, at 9 kpc (6 K-band disk scale lengths), 
M/Ly ~ 1 .5 (see Fig. 4). Since this increase runs counter 
to inside-out scenarios of galaxy formation, we must 
look for possible errors in our analysis. 

POSSIBLE PROBLEMS 

There are several possible sources of error in our anal- 
ysis. Some H II regions may be contaminating the PN 
sample and artificially lowering the PN velocity disper- 
sion. This is unlikely, since our photometric procedures 
should have eliminated virtually all such objects, but it 
must be considered. Similarly, though our analysis used 
an assumed value of 175 pc for the PN scale height, a 
25% change in this number does not significantly alter 
the results. Finally, our measurement of the disk matter 
scale length could be affected by a systematic change in 
the galaxy's internal extinction, or a breakdown in the 
various stability arguments. Our analysis suggests that 
such a breakdown would cause us to underestimate the 
disk scale length, not overestimate it, but this possibil- 
ity cannot be excluded (see C04). Further observations 
of more face-on galaxies are needed to check these as- 
sumptions. 

CONCLUSIONS AND FUTURE WORK 

For at least two decades, astronomers have been studying 
dark matter halos by assuming that disk mass-to-light 
ratios are constant. Our results for the PN kinematics of 
M33 suggests that this assumption is not valid. However, 
because M33 is inclined 56° to the line-of-sight, there 
are a number of caveats associated with our conclusions. 
To remove these uncertainties, we have begun studies 



REFERENCES 

1. Ashman, K.M. 1992, Pub. A.S.P., 104, 1 109 

2. Binney, J., & Tremaine, S. 1987, Galactic Dynamics, 
(Princeton: Princeton Univ. Press) 

3. Bizyaev, D., & Mitronova, S. 2002, Astr. Ap., 389, 795 

4. Bottema, R. 1993, Astr. Ap., 275, 16 

5. Bottema, R., van der Kruit, PC, & Freeman, K.C. 1987, 
Astr. Ap., 178, 77 

6. Carlberg, R.G. 1987, Ap. J., 322, 59 

7. Ciardullo, R., Durrell, PR., Laychak, M.B., Herrmann, 
K.A., Moody, K., Jacoby, G.H., & Feldmeier, J.J. 2004, 
Ap. J., 614, 167 

8. Ciardullo, R., Feldmeier, J.J., Jacoby, G.H., de Naray, 
R.K., Laychak, M.B., & Durrell, PR. 2002, Ap. J., 577, 
31 

9. Combes, F. 2002, New Astronomy Rev. , 46, 755 

10. Faber, S.M., & Gallagher, J.S. 1979, ARA&A, 17, 135 

1 1 . Feldmeier, J.J., Ciardullo, R., & Jacoby, G.H. 1996, Ap. J. , 
461, L25 

12. Gerssen, J., Kuijken, K., & Merrifield, M.R. 1997, 
MNRAS, 288, 618 

13. Gerssen, J., Kuijken, K., & Merrifield, M.R. 2000, 
MNRAS, 317, 545 

14. Herrmann, K.A., Ciardullo, R., & Vinciguerra, M. 2005, 
Bull. A.A.S., 205, 138.12 (2005). 

15. Jenkins, A., & Binney, J. 1990, MNRAS, 245, 305 

16. Kent, S.M. 1986, A. J., 91, 1301 

17. Merritt, D., & Sellwood, J. A. 1994, Ap. J., 425, 567 

18. Morosov, A.G 1980, Soviet Astron., 24, 391 

19. Morosov, A.G. 1981a, Soviet Astron., 25, 19 

20. Morosov, A.G. 1981b, Soviet Astron., 25, 421 

21. Palunas, P., & Williams, T.B. 2000, A.J., 120, 2884 

22. Sofue, Y., Koda, J., Nakanishi, H, & Onodera, S. 2003, 
Pub. A.S.J. , 55, 59 

23. Tilanus, R.P.J., & Allen, R.J. 1993, Astr. Ap., 274, 707 

24. Toomre, A. 1964, Ap.]., 139, 1217 

25. Toomre, A. 1966, in Notes on the 1966 Summer Study 
Program in Geophysical Fluid Dynamics at the Woods 
Hole Oceanographic Institution, (Woods Hole: Woods 
Hole Oceanographic Institute), 1 1 1 

26. van der Kruit, PC, Jimenez- Vicente, J., Kregel, M., & 
Freeman, K.C. 2001, Astr. Ap., 379, 374 

27. van der Kruit, PC, & Searle, L. 1981, Astr. Ap., 95, 105 

28. Villumsen, J.V. 1985, Ap. J., 290, 75 



